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Abstract 

Azimuthal asymmetries in exclusive electroproduction of a real photon from an unpolarized 
deuterium target are measured with respect to beam helicity and charge. They appear in the 
distribution of these photons in the azimuthal angle cj> around the virtual-photon direction, rel- 
ative to the lepton scattering plane. The extracted asymmetries are attributed to either the 
deeply virtual Compton scattering process or its interference with the Bethe-Heitler process. 
They are compared with earlier results on the proton target. In the measured kinematic region, 
the beam-charge asymmetry amplitudes and the leading amplitudes of the beam-helicity asym- 
metries on an unpolarized deuteron target are compatible with the results from unpolarized 
protons. 



Key words: DIS, HERMES experiment, GPD, DVCS, deuteron, unpolarized deuterium target 
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Fig. 1. Leading order Feynman diagrams for (a) deeply virtual Compton scattering and (b) the 
Bethe-Heitler process. 



1. Introduction 

Lepton-nucleon scattering experiments have long been an important tool in the de- 
tailed study of nucleon structure [I] . Two complementary approaches have contributed 
the most to our understanding of the nucleon. Elastic lepton-nucleon scattering has been 
exploited to extract nucleon form factors, which reveal how the electromagnetic nucleon 
structure differs from that of a point-like spin-1/2 particle. In another approach, Par- 
ton Distribution Functions (PDFs) are extracted from Deeply Inelastic Scattering (DIS). 
They represent distributions in the longitudinal momentum fraction carried by quarks 
and gluons in a nucleon moving with "infinite" momentum. PDFs and form factors 
present only one-dimensional pictures of nucleon structure. In recent years, a more com- 
prehensive multi-dimensional description of the nucleon has emerged in the framework of 
Generalized Parton Distributions (GPDs) [21314] . Their dependence on three kinematic 
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Fig. 2. Definition of the azimuthal angle <j> between the lepton scattering and photon production planes. 
Note that the azimuthal angle defined in this work differs from that used in Ref. [15] : <j> = tt — tftrjgj . 



quantities in addition to their evolution with the hard scale of the process carries informa- 
tion on two-parton correlations and quark transverse spatial distributions |5l6l7l8l9ll0j . 
GPDs embody PDFs as limiting cases, while elastic form factors appear as certain GPD 
moments. Other moments are connected with the total parton angular momentum con- 
tribution to the nucleon spin via the Ji relation [4] . 

GPDs can be constrained by measurements of hard exclusive leptoproduction of a pho- 
ton or meson in 'elastic' processes that leave the target intact. In Deeply Virtual Compton 
Scattering (DVCS), a quark absorbs a hard virtual photon, emits an energetic real photon 
and joins the target remnant (see Fig. [T](a)). DVCS is presently the only experimentally 
feasible hard exclusive process for which the effects of next-to-lcading order |11|12|13] 
and next-to-leading twist |14|15|16] are under complete theoretical control [17] . 

The final state of the DVCS process cannot be experimentally distinguished from that 
of the Bcthe-Heitler (BH) process, i.e., radiative elastic scattering (sec Fig.[T](b)). Hence, 
the two processes can interfere. Exclusive leptoproduction on a nucleon or nuclear target 
A of a real photon with four-momentum q' is denoted by 

e(k)+A(p)^e(k')+A(p') +1 (q'), (1) 

where k {k') and p {p') are the four-momenta of the incoming (outgoing) lepton and 
target, respectively. Averaged over the kinematic acceptance of the HERMES experiment, 
the BH cross section is much larger than that of the DVCS process. However, the BH 
cross section has a much weaker Q 2 dependence than the evolution of the DVCS cross 
section [4], so that in the HERMES energy range they can become comparable near 
Q 2 = 1 GcV 2 , with -Q 2 = q 2 = (k- k') 2 . 

Even in kinematic conditions where the DVCS process makes only a small contribution 
to the photon production cross section, its interference with the BH process provides 
access to the DVCS amplitudes through measurements of cross section asymmetries with 
respect to the charge and helicity of the incident lepton and the polarization of the target. 
These asymmetries appear in the distribution of the real photons in the azimuthal angle 
4>, defined as the angle between the lepton scattering plane, i.e., the plane defined by the 
incoming and outgoing lepton direction and the photon production plane spanned by the 
virtual and real photons (see Fig. [2]) . Significant azimuthal beam- helicity asymmetries in 
hard elcctroproduction of photons on the proton were first reported in Refs. [18119) . Later, 
asymmetries with respect to longitudinal [20)21] and transverse [22] target polarization, 
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as well as beam charge [33] and, with greater precision, beam helicity [23:25 26 27], were 
also measured on the proton. 

Measurements of azimuthal asymmetries for DVCS on nuclear targets [28] were ad- 
vocated as a useful source of information about partonic behavior in nuclei and nuclear 
binding forces |29] . If the target nucleus remains in its ground state the process is called 
coherent, while it is called incoherent if the nucleus is broken up. The dcutcron is a 
spin-1 nucleus, with implications for DVCS observables for the coherent reactions, which 
contribute mainly at very small values of the momentum transfer to the target. The 
asymmetries from the incoherent process involve mainly hard exclusive electroproduc- 
tion of a photon on the proton. The neutron contribution to the yield is typically small 
due to the suppression of the BH amplitude on the neutron by the small elastic electric 
form factor at low and moderate values of the momentum transfer to the target. 

This paper reports the first observation of azimuthal asymmetries with respect to beam 
helicity and charge for exclusive electroproduction of a real photon from an unpolarizcd 
deuterium target (e ^ d — > e ± J X). The dependence of these asymmetries on the kine- 
matic conditions of the reaction is also presented and certain asymmetry amplitudes 
are compared with the corresponding amplitudes obtained on an unpolarized hydrogen 
target (e ± p -> e ± 7 X) at HERMES [27]. 

2. GPDs and DVCS 

2.1. Generalized Parton Distributions 

In the generalized Bjorkcn limit of large Q 2 at fixed values of the Bjorken scaling 
variable xb — Q 2 /{2p ■ q) and small squared four-momentum transfer t = (p — p') 2 to 
the target, the DVCS process can be described by the leading (handbag) diagrams in 
Fig. QJa). Here, the process factorizes |3ll2l30j into a hard photon-quark scattering part 
calculable in quantum electrodynamics, and a soft part describing the nucleon structure, 
which can be expressed in terms of GPDs |2|3|4] . 

Like PDFs, GPDs depend on x and on the factorization scale Q 2 . In addition, GPDs 
depend on a skewness variable £ and the Mandclstam variable t. The skewness £ repre- 
sents half the difference in the longitudinal momentum fractions of the quark before and 
after the scattering, while x is their mean value (following the convention of Ref. [4]). In 
leading order, £ is directly accessible as it is related to the Bjorkcn scaling variable xb 
by £ ~ xb/(2 — xb)- In contrast, x is not directly accessible in DVCS, and some observ- 
ables appear as cc-convolutions of GPDs. Hence x plays a role different from that of xb 
in inclusive DIS. GPDs evolve logarithmically with Q 2 in analogy with PDFs [2I3I4I31| . 
This dependence on Q 2 is omitted for simplicity in the following. 

DVCS on spin-1/2 targets, such as nucleons, is described by four leading-twist quark- 
chirality conserving GPDs for each quark flavour q (and also for the gluon g), namely the 
GPDs m, E q , H q and E q [15]. The GPDs H q and E q arc quark-helicity averaged whereas 
H q and E q are quark-helicity dependent. The GPDs H q and H q conserve nucleon- helicity 
while E q and E q are associated with a helicity flip of the nucleon. In contrast, the coherent 
process on spin-1 nuclei, such as the deuteron, requires nine GPDs [32] — Hf, if| , if|, 
H%, Iff, Hf, If |, if| and H\ — to describe all DVCS observables. In the forward limit of 
vanishing momentum difference between the initial and final hadronic state (t — > and 
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£ — ► 0), the GPD H q (x,0, 0) reduces to f?(x), the quark number density distribution, 
and H q (x,0,0) reduces to g\{x), the quark helicity distribution. Similarly, for spin-1 
targets the GPDs H\ , H\ and H$ reduce to the following parton densities in the forward 
limit: 

Hf(x,0,0) 
H!(x,0,0) 

where i x ) represents the number density of a {anti} quark with momentum fraction 

{x < 0} x > and positive [negative] helicity in a rapidly moving deutcron target with 
longitudinal spin projection A. The 'unpolarizcd' (polarization averaged) quark densi- 
ties q A are defined as q A (x) = q^(x) + q A _(x). While the probabilistic interpretation of 
polarization-averaged and polarization-difference structure functions fi(x) and gi(x) in 
terms of quark densities is similar to that in the spin-1/2 case, the tensor structure func- 
tion bi(x) does not exist for spin-1/2 targets. It has been measured in DIS on a polarized 
spin-1 target [33]. Both H3 and H§ are associated with the 5% D-wave component of 
the deuteron wave function in terms of nucleons |34j . H3 is related to isoscalar currents 
and probes the binding forces in the deuteron, and H§ involves a tensor term [32|35j , the 
analog of which has no relationship to any local current due to Lorentz invariance. 



_ q 1 {x) +q- 1 (x) +q a (x) 
3 

= ql(x)-q-\x) 

= q°(x)- ql{x)+q ~ 1{x) 



/?(*), (2) 
<??(*), (3) 
bi(x), (4) 



2.2. Deeply virtual Compton scattering amplitudes 

For a target of atomic mass number A, the cross section for the hard exclusive lepto- 
production of real photons is given by [35 36^ 

da _ x A e 6 |T| 2 

Ax A dQ 2 d\t\ d</> ~ 32 (2tt) 4 Q 4 ' [ ' 

where xa = Q 2 /{ZMav) is the nuclear Bjorken xb, where Ma is the target mass and 
v = p ■ q/MA, £ = ^xaMa/ \/~Q 2 ~, and |T| is the total reaction amplitude. 

As the final states of the DVCS and BH processes are indistinguishable, the cross 
section contains the square of the coherent sum of their amplitudes: 

\1~\ 2 = |7bh + 7bvcs| 2 = |7bh| 2 + |7bvcs| 2 + Tbvcs ^bh + ^dvcs ^bh ■ (6) 

s ' 
I 

Here, I denotes the BH-DVCS interference term. The BH amplitude is calculable to 
leading order in Quantum Electrodynamics (QED) using nuclear form factors measured 
in clastic scattering. 

The interference term I in Eq. [H] provides separate experimental access to the real and 
imaginary parts of the DVCS amplitude through measurements of various cross-section 
asymmetries as functions of the azimuthal angle (f> [36] . Each of the three terms of Eq. [6] 
can be written as a Fourier series in <f> [15] . which in the case of an unpolarizcd target 
reads 



G 



I^Hl^^pxgcBHcOsM, (7) 

|T DVCS | 2 = A DVCS x {c° vcs + J2 C ° VCS c ° s (^) + ^? VCS ^} , (8) 

n=l 

1 = ~ Vi^)V 2 {^) X { C " + 51 C ™ cos ( n ^) + A s n sin WO} • (9) 

Here, -ft^H, -^dvcSj and K\ are kinematic factors, eg denotes the lepton beam charge in 
units of the elementary charge, and A the helicity of the longitudinally polarized lepton 
beam. The squared BH and interference terms have an additional cos (f> dependence in the 
denominator due to the lepton propagators Vi((f>) and V2 {(f>) in the BH process j!5|36j . 

The Fourier coefficients c l n and s l n in Eq. [5] can be expressed as linear combinations of 
Compton Form Factors t) (CFFs) [35], which in turn are convolutions of the corre- 
sponding GPDs F q (x,t;,t) with the hard scattering coefficient functions [11112113] : 



(10) 



where the — {+} sign applies to F q = Hf, . . . , i/| ■ ■ • , -Hi } in the case of a spin-1 

target. The real and imaginary parts of the CFFs have different relationships to the flavor 
sum over the respective quark GPDs. To leading order in a s , 

3m{^(C,<)} = -^^ e 2(F«(e,C,i) T i ;l<? (-^^i)) • (11) 

<? 

Hence measurements of cross-section asymmetries with respect to the beam helicity di- 
rectly determine combinations of GPDs along the lines x = ±£. In contrast, the real 
parts of the CFFs involve the full interval in x and constrain the x dependence of GPDs 
through convolutions: 

1 1 



q 



pJ^x Fi(x,£,t) 



(12) 



to leading order in a s . Here, P denotes Cauchy's principal value. Since the x dependence 
of GPDs is thereby only weakly constrained, experimental asymmetries in beam charge 
must be compared to the predictions of various GPD models. 

At leading twist (twist-2), the coefficients c\ and s\ are related to the same combina- 
tion of GPDs. This is also true for the kinematically suppressed coefficient Cq cx — ^^-c\. 
The coefficients cf, and c\ are sensitive to the 'D-tcrm' 37 38 , which contributes only in 
the 'ERBL' region — £ < x < £ where quark GPDs have the characteristics of distribu- 
tion amplitudes for the creation of a quark-antiquark pair. It does not contribute in the 
complementary 'DGLAP' region \x\ > £, where quark GPDs describe the emission and 
reabsorption of an (anti-)quark in the infinite momentum frame, thereby having proper- 
ties analogous to the familiar (anti-)quark distribution functions. The D-term provides 
a convenient means of representing this profound difference in GPD properties between 
the two regions, while, e.g., the absorption of this contribution into the double distri- 
butions [213] would require the introduction of terms with unnatural divergence, having 
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a severity beyond representation by delta functions or their derivatives. In addition to 
c\, s\, and c\, the only other Fourier coefficient related to only twist-2 quark GPDs is 
c dvcs^ The coe gj cients c pvcs 5 s pvcs ; c i ; and s i appear at tne twist-3 level, while c!? vcs 

and c\ arise from the gluonic transversity operator |39|40|41] at twist-2 level. The high- 
est harmonics of the interference and squared DVCS terms may also receive a twist-4 
contribution [42] . 

2.3. Azimuthal cross section asymmetries 

The beam-hclicity asymmetries for a longitudinally (L) polarized lepton beam and 
an unpolarized (U) target, based on the difference and sum of yields for the two beam 
charges, respectively, are defined as 

A i , _ [da + ^) - da+-(0)] - [da— (0) - da— (0)] 
LuW ~ [dcr^M + <kr+-(0)] + [da— (cf,) + da— (cf>)} ' 1 ' 

.DVCS U} = [da+-(0) - da+-(0)] + [da— (0) - da— (0)] 
^ LU W " [da+- (</>) + da+^ (</>)] + [da— (0) + da— (</>)} 1 ( ' 

where — > (<— ) denotes positive (negative) beam hclicity and the superscript + (— ) corre- 
sponds to positron (electron) beam. These definitions serve to separate the sin(n</>) terms 
in Eqs.[8]and[9] Similarly, the beam-charge asymmetry (BCA) for an unpolarized beam 
scattering from this target is defined as 

^^ da+^+da-W (15) 

~~ [da+^(0) + da+^ (</))} + [da— (<£) + da— (</>)] ' 
In terms of the Fourier coefficients of Eqs. [7H9] these equations read as 

4(« = —if -A^ W) (16) 

En=o <£ H cos(n^) + A DVC s ELo c° vcs cos(n0) 

-4E^ S (0) = —if 2 ^vc^i VCSsh ^ , (17) 

yi W% W ELo ^ H cobM + A DV cs Eto ^ VCS cos(n0) 



_4 / .x = -Viwkw ELo c » cosM 

PT(fe) ELo <5 H cos(n^) + A DV cs ELo ^ VCS cos(n0) ' 

At leading twist (twist-2, twist-3, and twist-2, respectively in the preceding three equa- 
tions), and neglecting gluonic terms, they reduce to 

rs? sinr^ 



A i (d \ ^ — ; (19) 

Prfe)ELo^ H cos(n0) + A DV csc|? vcs ' 

^LU W " K RH BH pogf^ , ^' r DVCS ' W 
, / 6) ~ ViWhw (4 + ^008 0) 

PrfeyELo^cos^ + Apvcscr 03, 
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To the extent that the DVCS contributions to the common denominator can be ne- 
glected at HERMES kinematics, the lepton propagators Vi((j>) and Vi{^) cancel in 
Eqs. [T6l [TH [HI and [21] However, this approximation is not invoked in the following 
because it would be subject to substantial model uncertainty. 



2.4. From Compton form factors to asymmetries 



Measured asymmetries are used to constrain GPD models by direct comparison of the 
data with model predictions. However, it is instructive to consider certain approximations 
relating CFFs and thereby GPDs to observed asymmetries. (These approximations are 
not needed in the comparison of GPD model predictions with measured asymmetries.) 

For an unpolarized nucleon target, the photon- helicity-conserving amplitude is 
given at leading twist by a linear combination of the CFFs H, H and £, together with 
the Dirac and Pauli form factors F\ and F2 1151 : 



M 1 - 1 =F 1 H + 7 ^—{F 1 +F 2 )H-j^ T F 2 £, (22) 



where xn is the Bjorken variable for the nucleon and Mm is the nucleon mass. At small 
values of xn and —t, M. ' — F\ Ji for the proton. For the neutron, the term containing 
the CFF £ in Eq. [22] becomes substantial at large —t due to the relative magnitudes 
of the form factors F\ and F2 for the neutron. The_Jcading Fourier coefficients of the 
interference term can be approximated as s\ cx SaW 1 ' 1 and c\ oc SteM 1 ' 1 . To leading 
order in 1/Q and in HERMES kinematic conditions, 

4Lufa)«-^sin0, (23) 

A c (4>) oc ^^cos0. (24) 

For the coherent process on the deuteron, the relationship between the Fourier coeffi- 
cients and the GPDs is complicated. However, the coefficients can be expanded in powers 
of Xd, the Bjorken variable for the deuteron target, and r = t/(4M D ), where Mr> is the 
deuteron mass [35] ■ Then, to leading order in a s and 1/Q, ^lu(^) can ^ e expressed in 
terms of the imaginary part of the deuteron CFFs Hi, H3 and H5 and the deuteron 
elastic form factors [13] G\ and G3 (see Fig.[3]). The quantity |r| is typically about 0.003 
in the range of small —t where the coherent process is significant, extending up to values 
of r only as large as 0.01. However, as shown in Fig.[3J the magnitude of G3 exceeds that 
of G\ by more than one order of magnitude. Hence certain terms leading in r (but not 
Xd) are retained. Defining 

= 3GWi - 2t[G 1 H 3 + G 3 (H 1 - ±H 5 )} + Ar 2 G 3 H 3 



ZG\ -4rGiG 3 +4t 2 G: 



the kinematic expansion yields 

x D (2-y)^(l-y) 
2-2y + y 

where y = p • q/(p ■ k). 



Av(<t>) Z Q \, 2 Qta^ 1 sin,/), (26) 
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The relative contributions of those terms are also shown in Fig. 02 they are less than 
10% at —t < 0.03 GeV 2 . When these terms are neglected, Eg. [251 becomes 



AuW- 2-2y + ^ GT sln ^ ( 27 ) 

The deuteron Ac (4>) is related to the real part of the same linear combination of CFFs 
appearing in the deuteron A\jj (</)): 



x D J-Qj{l-y) 

— ReD^cos^ (28) 



~ — — cos0. (29) 

y d 

For the coherent process on the deuteron, the leading term in the expansion of coeffi- 
cients s\ and c\ lead respectively to Eqs. [27l and [29l which arc analogous to Eqs. [23l and 
l24l for scattering on the nucleon. 



10 



Table 1 

The beam charge and polarization as well as the integrated luminosity in pb — 1 of the data sets used for 
the extraction of the various asymmetries on the unpolarized deuterium target. 

Beam Beam Luminosity 

Year Charge Polarization [pb — 1 ] 







A = -1 A = +1 


A = -1 A 


= +1 


1996 


e+ 


0.516 




43.9 


1997 


e+ 


-0.511 


53.1 




1998 


e — 


-0.307 


24.1 




1999 


e+ 


-0.552 0.418 


0.9 


5.1 


2000 


e+ 


-0.584 0.552 


29.7 


9.0 


2005 


e~ 


-0.355 0.377 


66.3 


65.7 


Sum 






174.1 


123.7 



3. The HERMES experiment 

A detailed description of the HERMES experiment can be found in Ref. 02] . A longi- 
tudinally polarized positron or electron beam of 27.6 GeV energy was scattered from an 
unpolarized deuterium gas target internal to the HERA lepton storage ring at DESY. 
The lepton beam was transversely polarized via the asymmetry in the emission of syn- 
chrotron radiation (Sokolov-Ternov effect) [45] in the arcs of the HERA storage ring. The 
transverse beam polarization was transformed locally into longitudinal polarization by 
a pair of spin rotators located before and after the experiment [46] . The hclicity of the 
beam was typically reversed approximately every two months. 

The beam polarization was continuously monitored by two Compton backscattering 
polarimeters |471 ; 4 8 . The average values of the beam polarization for various running 
periods are given in Table Q] the average fractional systematic uncertainty was 2.4%. 
The scattered leptons and produced particles were detected in the polar angle range 
0.04 rad < 8 < 0.22 rad. The lepton trigger required a coincidence of signals from 
scintillator hodoscope planes and the local deposition of a minimum energy of 3.5 GeV 
in the electromagnetic calorimeter. Lepton identification was accomplished using the 
transition-radiation detector, the preshowcr scintillator counter, and the electromagnetic 
calorimeter. The average lepton identification efficiency was at least 98% with hadron 
contamination that was less than 1%. Photons were identified by the detection of energy 
deposited in the calorimeter and preshowcr counter with no associated chargcd-particlc 
track. 

4. Event selection and yield distributions 

The data sets used in the extraction of the various asymmetries reported here are given 
in Table [1] In this analysis, it was required that events contained exactly one charged- 
particle track consistent with being the scattered beam lepton, and a single cluster in 
the calorimeter with an energy deposit E 1 > 5.0 GeV and with no associated charged 
track. The following requirements were imposed on the event kinematics: 1 GeV 2 < 
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Fig. 4. The measured distribution (points) of electroproduced real-photon events versus the squared 
missing mass . The solid curve represents a Monte Carlo simulation including coherent and incoher- 
ent BH and DVCS processes, the BH processes with the excitation of resonant final states (represented 
separately by the dashed-dotted curve), and the semi-inclusive background (dashed curve). The simula- 
tions and data are both normalized to the number of DIS events. The region between the two vertical 
lines indicates the selected exclusive events. 

Q 2 < lOGeV 2 , W% > 9GeV 2 , v < 22GeV and 0.03 < x N < 0.35, where W% = 
M$j + 2M N v — Q 2 , x N = Q 2 /(2M N v), and v = p ■ q/M N . The nucleonic (proton) mass 
M;v was used in all kinematic constraints on event selection even at small values of 
— t, where coherent reactions on the deuteron are dominant, because the experiment 
did not distinguish between coherent and incoherent scattering and the latter dominates 
over most of the kinematic range. Monte Carlo studies have shown that this choice has 
little effect on the extracted asymmetries [49]. In order to reduce background from the 
decay of neutral mesons, the angle between the laboratory 3-momcnta of the real and 
virtual photons was limited to # 7 » 7 < 45mrad. The minimum angle requirement # 7 * 7 > 
5 mrad was chosen according to Monte Carlo studies to be compatible with the effects of 
instrumental resolution in determination of <f>. 

'Exclusive' single-photon events were selected by requiring the squared missing mass 
M\ to be close to the squared nucleoli mass Mjy, where M\ is defined as M\ = 
(q + Pjv — <?') 2 w hh P]y = (M/v, 0, 0, 0). Due to the finite resolution of the spectrometer 
and the calorimeter, M|- may be negative. In Fig. 01 the squared missing mass distribu- 
tion of the selected events is compared with the predictions of Monte Carlo simulations 
of processes that contribute to both signal and background. One of the simulations uses 
an exclusive-photon generator for the BH and DVCS processes, including coherent and 
incoherent reactions as well as the excitation of resonant final states (a category known 
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0.2 0.4 0.6 

-t [GeV 2 ] 

Fig. 5. Distribution in —t of events selected in the exclusive region of M^- The points represent ex- 
perimental data while the continuous curve represents the simulation of real-photon production for all 
exclusive final states including resonances. Background from it decay is not included. The dotted and 
dashed curves represent the BH plus DVCS contributions of the coherent and incoherent elastic process, 
respectively. The dash-dotted curve shows the resonant BH contributions. The simulations and data arc 
both normalized to the number of DIS events. 

as associated production). The DVCS simulation for incoherent reactions on the proton 
is based on Ref. [50], while that for coherent reactions on the deuteron is based on the 
model from Ref. [35] . Most of the background in the vicinity of the exclusive peak comes 
from the decay of neutral pions. The dominant source of neutral pions is semi-inclusive 
DIS, 7*7V — > ir°X — > jjX, which is simulated using the Lepto event generator [5T] 
with a set of Jetset [52] fragmentation parameters tuned for HERMES kinematic con- 
ditions [53] , In this simulation, the photon originates mainly from decay of 7r°s from DIS 
fragmentation. Incoherent exclusive ir a production, 7*iV — > 7r°iV, was simulated using 
an exclusive Monte Carlo event generator based on the GPD models of Ref. [54] and 
was found to be negligible [J9 55j. HERMES data support this estimate [53]. The Monte 
Carlo yield exceeds the data by approximately 2% in the exclusive region. This may be 
due to the contribution of the DVCS process in the simulation of both coherent and 
incoherent processes, which is highly model-dependent and can vary between 10% and 
25% [55] for the incoherent processes. On the other hand, radiative effects not included 
in the simulation would move events from the peak to the continuum [57] , 

Events were selected in the 'exclusive region ', defined as -(1.5) 2 GeV 2 < M\ < 
(1.7) 2 GeV 2 to minimize background from DIS fragmentation while maintaining rea- 
sonable efficiency [58] . 
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As the recoiling target nucleon or nucleus was undetected, the Mandelstam variable 
t must be reconstructed from the measured kinematics of the scattered lepton and the 
detected photon. The resolution in the photon energy from the calorimeter is inadequate 
for a precise determination of t. Hence for events selected in the exclusive region in M^, 
the final state is assumed to be exclusive, leaving the target intact, thereby allowing t 
to be reconstructed with improved resolution using only the photon direction and the 
lepton kinematics [2"3"] : 

t = -Q 2 -2v{v- y/^ 2 + Q 2 cos 6> 7 * 7 ) ^ 

1 + JJ7 _ ^ + Q2 C0S6I 7*7 ) 

The further restriction —t < 0.7 GeV 2 is imposed in the selection of exclusive events in 
order to reduce background from the decay of neutral mesons. 

The t distribution of events for the deuterium target is shown in Fig. [5] and compared 
with the Monte Carlo simulations discussed above. The simulated contributions of co- 
herent and incoherent processes on the deuteron are also shown separately. Coherent 
scattering on the deuteron occurs preferentially at small values of —t. The Monte Carlo 
simulation shows that requiring -t < 0.06 GeV 2 enhances the mean fractional contribu- 
tion of the coherent process from 20% to 40% in the HERMES spectrometer acceptance. 
Requiring —t < 0.01 GeV 2 can further enhance the coherent contribution to 66%, but 
only at the cost of a rapidly decreasing yield. In Sections 16.11 and 16.21 the first two —t 
bins covering the range 0.00 — 0.06 GeV 2 will provide a measure of coherent effects; in 
Section [6~3l an attempt is made to isolate the coherent contribution. 



5. Analysis of the data 



5.1. Extraction of azimuthal asymmetry amplitudes 



The distribution of the expectation value of the yield for scattering a polarized lepton 
beam from an unpolarized deuterium target is given by 

(M)(Pt, e t , 4>) = C {P t , e t ) ti(e t , <j>) <tuu(0) 

x [1 + FMEu CS (^) + eeAc(<t>) + e*PMLuM] ■ (31) 

Here, C denotes the integrated luminosity Pi the longitudinal beam polarization, rj the 
detection efficiency, and auu (</>) the cross section for an unpolarized target averaged over 
both beam charges and both beam hclicitics, which can be expressed as 

fA\ XD 1 

gUuW = 32W VI + e* 

The asymmetries A\jj (0) , „4lij CS (4>) > an( i (4>) are related to the Fourier coefficients 
appearing in Eqs. [THll as illustrated by Eqs. [T6HT81 In analogy to the expansion of 
the cross section in Eq. 7-9, these asymmetries arc also expanded in terms of the same 
harmonics in <f>: 
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Al v (</>) * £ A*#*> sin(n0) + , (33) 



n=l 

^CS (0) „ ^in Dvcg gin ^ + ^cosjO^ ^ (34) 
3 

^ c (</))~^A™ sM) cos(n0), (35) 

where the approximation is due to the truncation of the in general infinite Fourier series 
caused by the azimuthal dependences in the denominators of Eqs. I16H18I 

For each kinematic bin in —t, xb, or Q 2 , the sets of azimuthal asymmetry amplitudes 
^LUI > ^lu'dvCS ano - A^ B ^ n ^\ hereafter called 'asymmetry amplitudes', are simulta- 
neously extracted from the observed exclusive sample using the method of maximum 
likelihood (described in detail in Ref. [H]). Although these asymmetry amplitudes differ 
somewhat from the coefficients given in Eqs. [7HH] and Eqs. [T6HT8I they are well defined 
and can be computed in various GPD models for direct comparison with the data. Note 
that in Eqs. 1331 and 1341 an additional constant term (n = 0) was introduced as a con- 
sistency test. These terms must vanish as they are parity violating. Removing these 
constant terms or also introducing additional harmonic terms in the fitting procedure do 
not influence results for other asymmetry amplitudes [59) . 

5.2. Background corrections and systematic uncertainties 

In each kinematic bin, the results from the maximum likelihood fit are corrected for 
photon background arising from semi-inclusive production of neutral mesons, mainly 
pions. A corrected asymmetry amplitude is obtained as 

. _ A'aw /semi ' ^scmi 

-^corr — Z 7 ■ (dO) 

-L /scmi 

Here, A ram stands for the extracted raw asymmetry amplitude, and / sem i and A sem i the 
fractional contribution and corresponding asymmetry amplitude of the semi-inclusive 
background. This fraction is obtained from a Monte Carlo simulation (see Section[4|) and 
ranges from 1% to 11%, depending on the kinematic conditions. As the semi- inclusive 
process is only very weakly beam-charge dependent, its asymmetry with respect to the 
beam charge or to the product of the beam charge and the beam polarization is assumed 
to be zero. The asymmetry of the semi-inclusive 7r° background with respect to only 
the longitudinal beam polarization is extracted from experimental data by requiring two 
photons to be detected in the calorimeter with an invariant mass between 0.10 GeV 
and 0.17 GeV and with no associated charged tracks. The restriction on the energy 
deposition in the calorimeter of the less energetic cluster is relaxed to 1 GeV to improve 
the statistical precision. The fractional energy z = E^/v of the reconstructed neutral 
pions is required to be larger than 0.8. After applying the correction of Eq. [36l the 
resulting asymmetry amplitudes are expected to originate from elastic (coherent), and 
incoherent photon production possibly including nuclcon excitation. 

The combined contribution to the systematic uncertainty from detector acceptance, 
smearing, finite bin width, and alignment of the detector elements with respect to the 
beam is determined from a Monte Carlo simulation using the GPD model described 
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Table 2 

The main contributions to the systematic uncertainty of extracted asymmetry amplitudes of interest, 
averaged over the full kinematic range. Not included is a 2.4% scale uncertainty of the beam-hclicity 
asymmetries due to the beam polarization measurement. 



Amplitude 


M\ shift 


Background corr. 


Acceptance, smearing, bin width, alignment 


.cos(O0) 


0.001 


0.001 


0.014 


A COS (ft 

A c 


< 0.001 


0.002 


0.023 


,sin</> 
A LU,I 


< 0.001 


0.004 


0.031 


^LU.DVCS 


0.002 


0.006 


0.003 



in Ref. [60]. Note that a mistake has been found in this GPD model [61]; however, the 
model described previously reported HERMES beam-charge [22] and preliminary (single- 
charge) beam-helicity asymmetries well [62] and thus is considered to be adequate for 
systematic studies. In each bin, the systematic uncertainty is taken as the difference be- 
tween the model prediction at the mean kinematic value of that bin and the respective 
amplitude extracted from the reconstructed Monte Carlo data. The dominant contribu- 
tions to the total systematic uncertainty are those from the detector acceptance and finite 
bin width. Further sources of uncertainty are associated with the background correction 
and a relative shift of the M\ spectra between the data samples from various running 
periods [22]. The contributions to the systematic uncertainty are added in quadrature. 
The main contributions for asymmetry amplitudes of interest are given in Table [21 Not 
included is any contribution due to additional QED vertices, as the most significant of 
these has been estimated to be negligible [6"5] . 

6. Results 

6.1. Results on beam- charge and beam-helicity asymmetries for an unpolarized 
deuterium target 

The asymmetry amplitudes are shown in Figs.[6H8]as a function of —t, xn, or Q 2 . While 
the variable xd would be the appropriate choice when presenting experimental results 
for pure coherent scattering, the nucleonic Bjorken variable xn is the practical choice 
in this case where incoherent scattering dominates over most of the kinematic range. 
The variables xn and Q 2 are strongly correlated due to the experimental acceptance. 
The 'overall' results in the left columns correspond to the entire HERMES kinematic 
acceptance. Figure [6] shows the amplitudes A£j , which are related to beam charge 
only, and Fig. [7] shows the amplitude -^lu'dvcs' which is related to beam helicity only, 
and the amplitudes A^u i \ which are related to both. All amplitudes are listed in 
Table [3] with the mean kinematic values of each bir l 17 1 . 

Of special interest is the asymmetry amplitude A^ s ^ , which is sensitive to the GPD Hi 
(H) for the coherent (incoherent) process in HERMES kinematic conditions (see Eqs.l29l 
and I24j) . The present data indicates that this amplitude increases with increasing —t. 



17 These results for only four bins in — t, XjVj or Q 2 , i.e., a binning used in previous HERMES pa- 
pers 1221231 . are available in the Durham database. 
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overall -t [GeV 2 ] x N Q 2 [GeV 2 ] 

Fig. 6. Amplitudes of the beam-charge asymmetry, which are sensitive to the interference term, in bins 
of — t, Xff, or Q 2 . The squares represent the results from the present work. The error bars (bands) 
represent the statistical (systematic) uncertainties. The finely (coarsely) hatched bands are theoretical 
calculations for incoherently combined proton and neutron targets, using variants of a double-distribution 
model 1 541661641 with the VGG Regge (VGG Factorized) ansatz for GPDs. The lowest panel shows the 
simulated fractions of coherent and resonant production. 

The amplitude ^4™ s (°^ in Fig. [SJ which is expected to relate to the same combination of 
GPDs as does A™ s ^, shows similar behaviour but with opposite sign, as expected [15]. 
The other two amplitudes A C q S ^ 2 ^ and A C q S ^\ related to twist-3 GPDs and the gluon 
transversity operator, respectively (see Section [2. 2[) . are consistent with zero. 

The fractional contributions to the yield from the coherent processes and from processes 
with excitation of resonant final states are presented in the bottom row of Fig. [6] (see 
also Table HJ, as obtained from the Monte Carlo simulation using the exclusive- photon 
generator mentioned in Section 3J Note that these fractional contributions are subject to 
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10 2 10" 1 10" 1 1 10 

overall -t [GeV 2 ] x N Q 2 [GeV 2 ] 

Fig. 7. The first row shows the sin <f> amplitude of the beam-helicity asymmetry that is sensitive to 
the squared DVCS term, in bins of — t, xjv, or Q 2 . Correspondingly, the second (third) row shows the 
sin (f> (sin 20) amplitude of the beam-helicity asymmetry sensitive to the interference term. All symbols 
are defined as in Fig. El There is an overall 2.4% scale uncertainty arising from the uncertainty in the 
measurement of the beam polarization. 



considerable model dependence. 

Figure [7] shows amplitudes of beam-helicity asymmetries, with the charge-averaged 
case related to the squared DVCS term in the upper row and the charge-difference case 
related to the interference term in the other rows. The amplitude ^lu'dvcS' wmcn is 
related to twist-3 GPDs, is found to be consistent with zero. Like the amplitude A^ s , the 
amplitude A^y j is also sensitive to the GPD Hi [H] for the coherent [incoherent] process, 
although these two asymmetries reveal different aspects of the (real) GPD, selected by 
different convolutions with (complex) hard scattering amplitudes. While the amplitude 
A^ s * is related to the real part of the CFF TCi [H] , the ^Lut amplitude is proportional to 

the imaginary part and shows significant negative values. The amplitude appears 
at twist-3 level, but nevertheless it shows a value which is non-zero and positive by 1.7 
standard deviations of the total experimental uncertainty. Figure [8] shows the amplitudes 
that are forbidden by parity conservation but were included in the fit as a consistency 
test. They are consistent with zero. 

The two hatched bands in Figs. [6] and [7] are theoretical calculations for the incoherent 
process, based on two different ansatze for modeling GPDs [BJ in the VGG model [53] 
(the coherent process will be considered in Section 6.3.). In this model, a GPD is written 
as a double distribution |2|3| complemented by a D-term [37I38| : 

- In the 'factorized ansatz' (VGG Fact.), the dependences on t and (x,£) are uncor- 
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Table 3 

Results for azimuthal Fourier amplitudes of the asymmetries with respect to the beam charge and helicity 
for the exclusive sample. 



kinematic bin {— t) 




.cos (0*) 

A c 


-cos 




.cos (24,) 


■COB (30) 

A C 


[GoV 2 ] 


[GeV 2 ] 


±#stat ± <5 sys t 


±<5 s tat 


± ^syst 


±5 sta t ± 5 syB t 


±<^stat ± <5 sys t 


overall 0.13 


0.10 2.5 


-0.028 ± 0.010 ± 0.014 


0.067 ± 0.015 


± 0.023 


-0.007 ± 0.014 ± 0.016 


0.005 ± 0.014 ± 0.001 



0.00-0.03 


0.02 


0.07 


1.7 


0.03-0.06 


0.04 


0.09 


2.2 


0.06-0.10 


0.08 


0.10 


2.-1 


0.10-0.20 


0.14 


0.11 


2.7 


0.20-0.35 


0.26 


0.12 


3.1 


0.35-0.70 


0.46 


0.11 


3.5 



0.03-0.06 


0.12 


0.05 


1.3 


0.06-0.08 


0.10 


0.07 


1.8 


0.08-0.10 


0.11 


0.09 


2.3 


0.10-0.13 


0.13 


0.11 


2.0 


0.13-0.20 


0.17 


0.16 


4.(1 


0.20-0.35 


0.23 


0.24 


(i.l 



-0.004 ± 0.023 ± 0.003 0.056 ± 0.035 ± 0.011 0.003 ± 0.032 ± 0.008 -0.007 ± 0.032 ± 0.004 
-0.020 ± 0.024 ± 0.014 0.034 ± 0.033 ± 0.007 -0.013 ± 0.033 ± 0.001 -0.041 ± 0.034 ± 0.005 
-0.008 ± 0.024 ± 0.020 0.045 ± 0.034 ± 0.023 0.014 ± 0.035 ± 0.017 0.059 ± 0.035 ± 0.010 
-0.029 ± 0.021 ± 0.024 0.085 ± 0.030 ± 0.027 -0.026 ± 0.029 ± 0.023 -0.006 ± 0.029 ± 0.006 
-0.067 ± 0.028 ± 0.018 0.093 ± 0.039 ± 0.022 -0.006 ± 0.038 ± 0.050 0.044 ± 0.037 ± 0.002 
-0.066 ± 0.042 ± 0.029 0.114 ± 0.064 ± 0.057 -0.015 ± 0.056 ± 0.049 -0.007 ± 0.055 ± 0.008 

-0.052 ± 0.026 ± 0.003 0.092 ± 0.040 ± 0.036 0.024 ± 0.031 ± 0.002 0.024 ± 0.030 ± 0.005 
-0.024 ± 0.022 ± 0.017 0.049 ± 0.031 ± 0.028 0.004 ± 0.029 ± 0.014 -0.027 ± 0.030 ± 0.001 
-0.030 ± 0.025 ± 0.016 0.056 ± 0.036 ± 0.023 -0.014 ± 0.035 ± 0.013 -0.008 ± 0.035 ± 0.008 
0.011 ± 0.026 ± 0.030 0.039 ± 0.037 ± 0.030 0.004 ± 0.037 ± 0.029 0.068 ± 0.037 ± 0.004 
-0.021 ± 0.028 ± 0.007 0.070 ± 0.040 ± 0.031 -0.051 ± 0.040 ± 0.011 0.000 ± 0.038 ± 0.002 
-0.013 ± 0.052 ± 0.055 0.136 ± 0.074 ± 0.022 -0.091 ± 0.069 ± 0.039 0.000 ± 0.069 ± 0.008 



1.0 


■1.4 


0.09 


0.05 


1 


2 


-0.032 ± 0.022 ± 0.016 


0.077 ± 0.033 ± 0.037 


0.000 ± 0.029 ± 0.007 


-0.030 ± 0.030 ± 0.004 


1.4 


■1.8 


0.10 


0.07 


1 


(i 


-0.050 ± 0.026 ± 0.017 


0.100 ± 0.037 ± 0.016 


0.000 ± 0.034 ± 0.018 


0.021 ± 0.034 ± 0.006 


1.8 


2.4 


0.12 


0.09 


2 


1 


-0.031 ± 0.023 ± 0.017 


0.025 ± 0.032 ± 0.035 


-0.035 ± 0.033 ± 0.015 


0.074 ± 0.034 ± 0.009 


2.4-3.2 


0.14 


0.11 


2 


8 


-0.021 ± 0.024 ± 0.026 


0.106 ± 0.038 ± 0.013 


0.045 ± 0.036 ± 0.021 


-0.042 ± 0.035 ± 0.003 


3.2 


■4.5 


0.16 


0.14 


3 


8 


-0.010 ± 0.027 ± 0.014 


0.026 ± 0.037 ± 0.029 


0.018 ± 0.037 ± 0.015 


-0.003 ± 0.037 ± 0.004 


4.5- 


10.0 


0.23 


0.20 


5 


8 


-0.010 ± 0.032 ± 0.035 


0.055 ± 0.046 ± 0.023 


-0.095 ± 0.046 ± 0.025 


0.013 ± 0.045 ± 0.003 



kinematic bin ( — t) (^jv) {Q ) 
[GeV 2 ] [GeV 2 ] 



LU.DVCS 
±£stat ± <5 sys t 



LU,I 

±5stat ± ^syst 



.sin (2<£) 
LU,I 
±<5stat ± 6 sys t 



0.13 0.10 2.5 



-0.007 ± 0.033 ± 0.007 -0.192 ± 0.035 ± 0.031 0.073 ± 0.031 ± 0.012 



0.00-0.03 


0.02 


0.07 


1.7 


0.03-0.06 


0.04 


0.09 


2.2 


0.06-0.10 


0.08 


0.10 


2.4 


0.10-0.20 


0.14 


0.11 


2.7 


0.20-0.35 


0.26 


0.12 


3.1 


0.35-0.70 


0.46 


0.11 


3.5 



-0.042 ± 0.074 ± 0.011 -0.296 ± 0.104 ± 0.006 0.056 ± 0.071 ± 0.011 
-0.101 ± 0.077 ± 0.013 -0.167 ± 0.084 ± 0.008 0.034 ± 0.072 ± 0.009 
0.032 ± 0.080 ± 0.032 -0.064 ± 0.081 ± 0.010 0.114 ± 0.076 ± 0.032 
0.018 ± 0.068 ± 0.009 -0.215 ± 0.071 ± 0.016 -0.022 ± 0.065 ± 0.013 
0.095 ± 0.087 ± 0.009 -0.286 ± 0.095 ± 0.008 0.206 ± 0.085 ± 0.024 
-0.029 ± 0.118 ± 0.035 0.003 ± 0.122 ± 0.005 0.133 ± 0.124 ± 0.030 





0.03-0.06 


0.12 


0.05 


1.3 




0.06-0.08 


0.10 


0.07 


1.8 




0.08-0.10 


0.11 


0.09 


2.3 


H 


0.10-0.13 


0.13 


0.11 


2.0 




0.13-0.20 


0.17 


0.16 


4.0 




0.20-0.35 


0.23 


0.24 


(i.l 



-0.007 ± 0.064 ± 0.021 -0.197 ± 0.083 ± 0.061 0.080 ± 0.066 ± 0.015 
0.012 ± 0.069 ± 0.018 -0.286 ± 0.096 ± 0.032 0.084 ± 0.067 ± 0.014 
0.041 ± 0.080 ± 0.025 -0.017 ± 0.080 ± 0.031 -0.018 ± 0.075 ± 0.010 
-0.056 ± 0.084 ± 0.033 -0.212 ± 0.090 ± 0.023 0.060 ± 0.080 ± 0.013 
-0.109 ± 0.090 ± 0.037 -0.189 ± 0.093 ± 0.020 0.029 ± 0.083 ± 0.002 
0.222 ± 0.160 ± 0.053 -0.313 ± 0.161 ± 0.032 0.444 ± 0.163 ± 0.032 



1.0-1.4 


0.09 


0.05 


1 


2 


-0.028 ± 0.068 ± 0.035 


-0.208 ± 0.082 ± 0.060 


0.052 ± 0.065 ± 0.011 


1.4-1.8 


0.10 


0.07 


1 


6 


0.175 ± 0.079 ± 0.030 


-0.222 ± 0.087 ± 0.049 


0.127 ± 0.077 ± 0.019 


1.8-2.4 


0.12 


0.09 


2 


1 


-0.108 ± 0.076 ± 0.020 


-0.124 ± 0.077 ± 0.029 


-0.011 ± 0.071 ± 0.007 


2.4-3.2 


0.14 


0.11 


2 


8 


0.005 ± 0.083 ± 0.023 


-0.244 ± 0.091 ± 0.037 


0.054 ± 0.077 ± 0.010 


3.2-4.5 


0.16 


0.14 


3 


8 


-0.045 ± 0.086 ± 0.037 


-0.169 ± 0.088 ± 0.022 


0.119 ± 0.082 ± 0.005 


4.5-10.0 


0.23 


0.20 


5 


8 


-0.038 ± 0.104 ± 0.010 


-0.233 ± 0.105 ± 0.006 


0.166 ± 0.100 ± 0.006 
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2.4 % scale uncertainty 





10" 

overall 
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Fig. 8. The cos(0</>) amplitudes (constant terms) that are included as a consistency test in the fit in 
Eqs, 1341 and 1331 All symbols arc defined as in Fig. [6] There is an overall 2.4% scale uncertainty arising 
from the uncertainty in the measurement of the beam polarization. 



related. The t dependence is written in accordance with proton elastic form factors. 
The (a;, £) dependence is based on double distributions [2] constructed from ordinary 
PDFs complemented with a profile function that characterizes the strength of the £ 
dependence; in the limit b — > oo of the profile parameter b, the GPD is independent of 
£ [66| . Note that b is a free parameter to be experimentally determined independently 
for valence and sea quarks. 
- The 'Regge ansatz' (VGG Rcgge) implements entanglement of the t dependence of 
the GPD with its dependence on x and £. This feature is inspired by the traditional 
interpretation of measurements of elastic diffractive processes in terms of Regge phe- 
nomenology [64], and finds further support in more recent phenomcnological con- 
siderations [67168] . This ansatz for GPDs hence uses for the t dependence of the 
double distributions a soft Regge-type parameterization oc |£|- Q (°)+ Q 1*1 with a' = 
0.8 GeV" 2 ... 0.9 GeV" 2 for quarks. 

Both theoretical calculations are averaged at the cross section level over incoherent 
processes on the proton and neutron in each kinematic bin. In both calculations the 
D-term is assigned the value zero. Earlier, it was found that inclusion of a D-term with 
any significant magnitude in the double-distribution model of Rcf. [65j employing several 
variants of Regge or factorized ansatze with any choice of profile parameters fails to 
describe the BCA amplitudes measured at HERMES on a hydrogen target [22I23| . The 
theoretical bands in Figs. [6] and [7] correspond to the range of values of the asymmetry 
amplitudes obtained by varying the profile parameters 6 va i and 6 sea between unity and 
infinity. The theoretical calculations based on the factorized ansatz fail to describe the 
t dependence of A™^ ^ and A^ s * as seen in Fig. [6] The calculations based on the 
Regge ansatz for GPDs are in good agreement with the t dependence of the measured 
asymmetry amplitudes with respect to the beam charge j\^ s ^\ Both ansatze predict 
that A™ s ^ decreases with increasing xjv, which is not seen in the data. Both ansatze 
undershoot the asymmetry amplitudes with respect to the beam helicity i ■ 
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Table 4 

Simulated fractional contributions for coherent and resonant processes in each kinematic bin. 



kinematic bin coherent resonant 





overall 


0.176 


0.174 




0.00 


- 0.03 


0.481 


0.064 




0.03 


- 0.06 


0.256 


0.110 


CN 

> 
1 


0.06 
0.10 


- 0.10 

- 0.20 


0.130 
0.053 


0.150 
0.206 




0.20 


- 0.35 


0.017 


0.289 




0.35 


- 0.70 


0.005 


0.387 




0.03 


- 0.06 


0.258 


0.161 




0.06 


- 0.08 


0.214 


0.160 




0.08 


- 0.10 


0.176 


0.173 




0.10 


- 0.13 


0.127 


0.184 




0.13 


- 0.20 


0.078 


0.203 




0.20 


- 0.35 


0.032 


0.198 




1.0 


- 1.4 


0.253 


0.133 


> 
o 


1.4 


- 1.8 


0.209 


0.154 


1.8 
2.4 


- 2.4 

- 3.2 


0.172 
0.150 


0.172 
0.193 


3.2 


- 4.5 


0.109 


0.219 




4.5 


10.0 


0.055 


0.237 



6.2. Comparison of the deuteron results with the HERMES results on beam-charge and 
beam-helicity asymmetries on the proton 

In Figs. [9r fTT1 the overall asymmetry amplitudes as well as their — t, xn, and Q 2 depen- 
dences, measured for the unpolarized deuterium target, are compared with the analogous 
results obtained from HERMES data on the proton [27] . 

The deuteron data include the coherent process e d — > e ± d r y, and the incoherent 
process e ± d — > e ± p n 7, where a nucleon may be excited to a resonance. The proton data 
include only e ± p — » e ± p 7 and the case with resonance excitation. Any difference that 
appears at small values of —t may be due to the coherent process. Monte Carlo simulations 
indicate that the incoherent process dominates for 0.06 GeV 2 < — t < 0.7 GcV 2 (sec 
Fig. [5]) • As shown in Figs.lMTTl the deuteron and proton results are found to be consistent 
in most kinematic regions. A possible difference in the last two —t bins of the amplitude 
^cos0 ^ gee k c ^ Q ^ c contributions of the neutron and its resonances. 

The proton and deuteron results for the amplitude A^™ 1 integrated over the acceptance 
differ by 2.5 times the total experimental uncertainties. This possible discrepancy is most 
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Fig. 9. Fourier amplitudes of the beam-charge asymmetry, which are sensitive to the interference term, 
in bins of —t, iejv, or Q 2 , extracted from deuteron data (squares) and from proton data (triangles). The 
points for deuterium are slightly shifted along the x-axis for visibility. The error bars (bands) represent 
the statistical (systematic) uncertainties. The hatched band is for the deuterium target. 



evident at large —t and large xn (or Q 2 ). Such a discrepancy would have no obvious 
explanation. 
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Fig. 10. Fourier amplitudes of the beam-helicity asymmetry that are sensitive to the squared DVCS term, 
in bins of — t, xjv, or Q 2 , extracted from deuteron data (squares) and from proton data (triangles). The 
error bars (bands) represent the statistical (systematic) uncertainties, which include all sources apart 
from the 2.4% (2.8%) scale uncertainty for the deuteron (proton) data due to the beam polarization. 
The hatched band is for the deuterium target. 

6.3. Estimates of the asymmetries from coherent scattering 



Estimates of the asymmetries for coherent scattering in the range —t < 0.06 GeV 2 , 
corresponding to the first two bins, were derived by correcting for the incoherent contri- 
butions of the proton and its resonances using the simulated fractional coherent contri- 
butions from Table [4l under the assumption that the asymmetries for these contributions 
are the same as those on the free proton. The simulated contribution of approximately 



7% from the process e : 



n 7 is estimated to have an effect on the asymmetries of 



less than 0.01. The extracted coherent asymmetries ^c°coh an d ^lu i con are found to be 



0.11 ± 0.07 (stat.) ± 0.03 (syst 
the average kinematic valuei 
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and -0.29 ± 0.18 (stat.) ± 0.03 (syst.), respectively, at 
(-*) = 0.03 GeV 2 , (x D ) = 0.04, and (Q 2 ) = 1.9 GeV 2 . 
These results for the coherent asymmetries are compared in Table [5] with model esti- 
mates using the models A, B, Bq, B, B' , and C of Rcfs. [15|35j . the main parameters 
of which are listed in Table [6] The model estimates are based on the double distribution 
ansatz |66j for nucleonic GPDs, combined with a factorized t dependence, and with the 
D-term set to zero. The nucleonic GPDs are combined using the impulse approxima- 
tion. The contribution of sea quarks is neglected in model Bq, while it is enhanced in 
model C by a choice of a smaller value of the parameter & sca , which increases the abso- 
lute value of the beam-helicity asymmetry amplitude i co h compared to model A. In 

model B' (J5), the GPD H% (i/5) is taken into account by arbitrarily equating it with 
Hi {Hi{x) — Hi(—x)). All other GPDs are kinematically suppressed and are set to zero. 
The models Bq and C were previously ruled out by the beam-helicity and beam-charge 
asymmetry measurements on the hydrogen target [18 19 22 2 3127) . 



Nucleonic Bjorken mjv is experimentally irrelevant for coherent scattering. 
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Fig. 11. Fourier amplitudes of the beam-helicity asymmetry that are sensitive to the interference term, 
in bins of —t, xjv, or Q 2 , extracted from deuteron data (squares) and from proton data (triangles). The 
error bars (bands) represent the statistical (systematic) uncertainties, which include all sources apart 
from the 2.4% (2.8%) scale uncertainty for the deuteron (proton) data due to the beam polarization. 
The hatched band is for the deuterium target. 



Table [5] also includes model predictions from Ref . [38] . This model is based on double 
distributions, where only the polarizations of the valence quarks are considered for the 
nucleonic GPDs. A factorized ansatz for the t dependence of the nucleonic GPDs is em- 
ployed and the strange quark contribution is neglected. Again the impulse approximation 
is used to combine the nucleonic GPDs, without including the particular contribution 
from the D-term. 

All models are consistent within two standard deviations in the total experimental 
uncertainty with the extracted results for ^Lut coh an< ^ ^c°coh' exce Pt for models Bo 
and that of Ref. [38], which disagree with the results of ^.^coh about 3.5 standard 
deviations. Here, it should be noted that predictions for the real part of the CFFs are 
subject to delicate cancellations |15j and hence are extremely sensitive to assumptions. 
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Table 5 

Experimental and theoretical values of the beam-helicity and beam-charge asymmetries for the coherent 
process on the deuteron. The theoretical predictions are for variants of the models of Ref. 115 1351 and a 
model from Ref. [38] . The experimental uncertainties do not account for the model dependence of the 
simulated fractional contributions of coherent and incoherent processes. 





Exp. value 








Model 










value ± <5 s tat ± (Ssyst 


A 


B 


Bo 


B 


B' 


C 


m 


* sin <ji 
LU,I,coh 

A COS (f> 

A C,coh 


-0.29±0.18±0.03 
0.11 ±0.07 ±0.03 


-0.44 
0.10 


-0.38 
0.09 


-0.16 
-0.17 


-0.37 
0.09 


-0.39 
0.09 


-0.58 
0.22 


-0.36 
-0.15 



Table 6 

Model parameter sets for the GPD H\ of the deuteron 15 35 . The t slope parameter B SEa is used mainly 
to change the normalization of the sea quark GPD Hi . 



deuteron Hi GPD 




Model 






Model parameters 


A 


B (B' , B) 


Bo 


C 


'<val 


1 


oo 


oo 


1 


^sca 


oo 


CO 




1 


Bsea [GeV- 2 ] 


20 


20 
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7. Summary 

Azimuthal asymmetries with respect to beam-helicity and beam-charge are measured 
for hard exclusive electroproduction of photons in deeply inelastic scattering off an unpo- 
larizcd deuterium target. The observed asymmetries are attributed to either the interfer- 
ence between the DVCS and the Bethe-Heitler processes or the pure DVCS process. The 
asymmetries are observed in the exclusive missing-mass domain — (1.5) 2 GeV 2 < M\ < 
(1.7) 2 GeV 2 . The dependences of these asymmetries on —t, Xn, or Q 2 are investigated. 
The results from the deuterium target include the coherent process e ± d — ► e d 7 and 
the incoherent process e ± d — > e pwy, where a nucleon may be excited to a resonance. 
For an unpolarized deuterium target, the leading Fourier amplitude of the beam-helicity 
asymmetry that is sensitive to the interference term is found to be substantial, but no 
significant t dependence is observed. The leading amplitude of the beam-charge asymme- 
try is substantial at large —t, but becomes small at small values of —t. The amplitudes of 
the beam-helicity asymmetry that are sensitive to the squared DVCS term are found to 
be consistent with zero. The data are able to discriminate among various GPD models. 

The measured asymmetry amplitudes from unpolarized deuteron and proton |27j tar- 
gets are consistent in most kinematic regions, except possibly for the leading amplitude 
of the beam-charge asymmetry in the last two —t bins, and the 'overall' value of A™^ 2 ^ . 

The beam-charge and beam-helicity asymmetry amplitudes for coherent scattering 
from the deuteron are extracted from the asymmetry amplitudes measured on unpolar- 
ized deuteron and proton targets. When compared to the GPD models of Refs. [15I35| . 
the results disfavor a large sea quark contribution while favoring a non-zero contribu- 
tion. The results disfavor the variants of the model of Refs. |15l35j that omit sea quark 
contributions, and also the model of Ref. [38 . 
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